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Determining hybrid content of heavy quarkonia using lattice 
nonrelativistic QCD 

Tommy Burch, Kostas Orginos*, and Doug Toussaint a 

a Department of Physics, University of Arizona, 
Tucson, AZ 85721 

Using lowest-order lattice NRQCD to create heavy meson propagators and applying the spin-dependent in- 
teraction, CBJm~^ ' ' a ^ var Y m g intermediate time slices, we compute the off-diagonal matrix element of the 
Hamiltonian for the quarkonium-hybrid two-state system. Diagonalizing this two-state Hamiltonian, the admix- 
ture of hybrid (\QQg)) in the ground state is found. We present results from a set of quenched lattices with an 
interpolation in quark mass to match the bottomonium spectrum. 



1. INTRODUCTION 

While experimental confirmation of the exis- 
tence of hybrid states - hadrons with a gluonic 
excitation (qqg, qqqg, etc.) - remains elusive, 
much has been done theoretically, both on and 
off the lattice, to determine possible consequences 
of their existence. In particular, mixing of (non- 
exotic) hybrids with normal hadronic states (with 
the same J PC ) is expected to occur, necessitating 
a redefinition of the true hadronic ground state in 
terms of its constituents: e.g., 



|T> = A s \bb) + A H \bbg) 



1) 



We offer the bottomonium example for a reason: 
it is this system which we study with our simula- 
tions. 

Since this system is composed of relatively 
heavy quarks, we use the nonrelativistic approx- 
imation to QCD (NRQCD) jl| when we perform 
the lattice simulations. NRQCD has been used 
before to study heavy hybrids [^J. The main dif- 
ference in the current study is the fact that we 
apply the lowest-order (in ^-) spin-dependent 

interaction [cgw^-a ■ B) "perturbatively," i.e., 
at a single time slice. Then, using different 
source and sink operators (quarkonium — > hy- 
brid, and vice versa), we are able to extract the 
off-diagonal matrix element of the Hamiltonian 
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for this quarkonium-hybrid two-state system. We 
then diagonalize this new Hamiltonian and find 
the admixture of hybrid within the true ground 
state |§. 

2. METHOD 

In our implementation of NRQCD, we use a 
time-step-symmetric form of the heavy quark evo- 
lution operator 0] : 



G{x,t + a) = 
aH 



1 - 



2n 



x (l-S t , t aSH)G(x,t), 



(2) 



with a value of n = 2 (more than sufficient for sta- 
bility, n > 2 m a) ' ^ or simplicity, we use only the 
lowest-order term of the heavy quark expansion 
in the diagonal part of our Hamiltonian: 



H = 



—D 1 
2m„ 



(3) 



where D is the covariant derivative. The spin- 
dependent term appears in the interaction 



SH = c B -^-5-B, 



(4) 



which is applied only at a single, intermediate 
time slice, t' . 

For our lattice mesons, we use an incoherent 
sum of point sources: at the source end, we 
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Table 1 

Meson operators. 



J 



Operator 



h S-wave 
1" S-wave 
++ P-wave 
1 ++ P-wave 

2++ P-wave Q{a l D 3 
0"+ Hybrid 

1 Hybrid 



QQ 
QcTiQ 
_ QaiDiQ 

Q £ ijk&jDkQ 
' - " a, 4 . 

Qa l B i Q 
QBiQ 



start with a given quark color and spin at all 
spatial points, without fixing the gauge; at the 
sink end, we sum over all the contributions where 
the quark and anti-quark are at the same spatial 
point. Since the lattices are not gauge-fixed, we 
expect the contributions from sources with the 
quark and anti-quark at different spatial points 
to average to zero. We combine the quark and 
anti-quark sources (propagators) with the appro- 
priate spin matrices and gauge fields to construct 
the meson operators at the source (and sink) time 
slices. The meson operators we use are displayed 
in Table @. 

Using identical source and sink operators, "un- 
mixed" meson (and hybrid) propagators are con- 
structed. We use the following form to fit the 
resulting correlators: 



C(t) = A e- m o' +A ie ' 



(5) 



For the "mixed" propagators, the source and sink 
operators differ (quarkonium — > hybrid, and vice 
versa). These are also fit with a two-exponential 
form. However, since the source and sink opera- 
tors arc different, and since there is the single in- 
teraction (SH) at the intermediate time slice (£') 
causing the mixing between these configurations, 
the first exponential in this fit can be expressed 



Cmix (^) 



A 1/2 A 1/2 



1H 



' 2m„ 



-a-B 



IS 



x e 



' i 0, S rc t 'g- m 0, Sn fc(*- t ') + 



(G) 



Knowing the amplitudes and masses from the un- 
mixed propagators and fitting this correlator in 



the region t > £', we extract the off-diagonal ma- 
trix element. We then repeat this calculation 
for larger and larger values of t' in search of a 
plateau, indicating the decay of excited-state con- 
tributions (e.g., from |2S) or |2H), depending on 
the source) to the mixed amplitude. 

3. RESULTS 

The heavy meson propagators are evaluated on 
a set of quenched, improved lattices (100 configu- 
rations with Symanzik 1-loop- improved gauge ac- 
tion, 20 3 x 64, (3 = 8.0; see Ref. {§). The lattice 
spacing is determined using the spin- averaged 1S- 
1P mass splitting for bottomonium; a -1 = (440 
MeV)/(aAM SP ) = 1590(30) MeV. We also con- 
struct non-zero- momentum 1 S-wave propaga- 
tors and use the resulting dispersion relation to 
determine the kinetic mass for this meson. This 
kinetic mass is calculated for two input values of 
the quark mass (m q a = 2.5 and 2.8) and an inter- 
polation is then performed to match this meson 
mass to the experimental value for the mass of 
the T (9.46 GeV). We thus arrive at a (lattice- 
regularized) physical value for the mass of the 
bottom quark (m b a = 2.71). 

Appearing in Fig. [J are the results for the off- 
diagonal matrix element, plotted as a function of 
the interaction time slice, t'. After diagonaliza- 
tion of the Hamiltonian, the hybrid content of 
the ground state is found and expressed in terms 
of a mixing angle (sin#). A plot of this quan- 
tity versus t' is shown in Fig. ||. As can be seen 
in both of these figures, a plateau is reached at 
rather low values of t' . Using the values at t' = 10 
(% 2 /d.o.f. < 1) and interpolating in the quark 
mass to the physical value, we find the following 
configurations for the T and rjb ground states: 

|T> = 0.99837(6) |1S(1~)> 

- 0.057(1) |1H(1— )), (7) 

1%) = 0.9953(1) |lS(0-+)> 

+ 0.097(2) |lH(0-+)>, (8) 

with the tree-level value c B = 1. 

In an attempt to determine c B nonperturba- 
tively, we perform additional runs with the inter- 
action, SH, applied at all intermediate time slices 
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Figure 1. Magnitude of the off-diagonal matrix 
element of the Hamiltonian (in lattice units) for 
the S-wave / Hybrid two-state system vs. the 
time slice, t', at which the interaction (8H) is 
applied (m q a = 2.5). 



and with the values of c B — 1 and 2. We then 
use the T — rj b mass difference (a quantity hav- 
ing a quadratic dependence upon the a ■ B term) 
to set the value of c 2 B . Unfortunately, as there is 
no experimental value for the mass of the rj b , we 
rely on potential model results p] which put this 
mass difference in the 30 — 60 MeV range. Com- 
bining this with our lattice results for the mass 
difference - AM T _ W = 18.4(4) and 71.3(2.1) 
MeV for c B = 1 and 2, respectively - suggests 
c B ~ 1.7 — 3.4. The resulting probability admix- 
tures of hybrids within the bottomonium ground 
states thus become 

|(1H|T}| 2 - 0.0054- 0.011, (9) 

|(lH|?7 fc )| 2 - 0.016 - 0.032. (10) 

Obviously, the uncertainty is dominated by the 
lack of a reliable value for the T — r\ b mass differ- 
ence. Note the factor of ~ 3 enhancement for the 
mixing in the channel, due to spin statistics 
jfjj. Interested readers can find a more detailed 
discussion of the results for the 1~~ channel in 
Ref. |. 



Figure 2. Mixing angle, sin(0), vs. the time slice, 
t' , at which the interaction term (SH) is applied 
(m q a = 2.5). 



The simulations are to be repeated on lattices 
with a different value of the coupling in order to 
attempt a continuum extrapolation. 
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